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Abstract 

We study the asymptotics of the spectral distribution for large empiri- 
cal covariance matrices composed of independent Multifractal Random Walk 
processes. The asymptotic is taken as the observation lag shrinks to 0. In 
this setting, we show that there exists a limiting spectral distribution whose 
Stieltjes transform is uniquely characterized by equations which we specify. 
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1. Introduction 

Since the seminal work of Marcenko and Pastur [8] in 1967, there has been growing 
interest in studying the asymptotics of large empirical covariance matrices. These 
studies have found applications in many fields of science: physics, telecommunica- 
tions, information theory and finance, etc... The main motivation of this work stems 
from finance: we refer to [9], [4] for a discussion on the applications of large empirical 
covariance matrices in finance and in particular in portfolio theory. In the context of 
finance, the study of the limiting spectral distribution of large empirical covariance 
matrices is of particular interest as it is a crucial statistical tool in identifying the 
different market modes (see again [9], [4] for a nice introduction to this topic). It is 
indeed very natural to try to identify common causes (or factors) that explain the 
dynamics of N quantities, which will be stock's prices in finance. We will denote by 
iV the number of stocks and by T the number of time intervals where we observed 
prices of the N stocks. In this setting, the Marcenko Pastur paper enables to deal 
with the case where stock prices follow independent Brownian motions. More pre- 
cisely, let us define a iV x T matrix Xn such that X^iij) is the realization of the 
return on the j-th time interval (of shrinking size 1/T) of stock number i by: 

X N (ij) = B, (£j - B t (1.1) 

where the Bi are i.i.d. standard Brownian motions. The empirical covariance matrix 
is now defined as Rn = X^X N . If Ai, . . . , Xn are the eigenvalues of Rn, the empirical 
spectral distribution of the matrix Rn is the probability measure defined by: 

1 - 
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The Marcenko and Pastur result states that, in the limit of large matrices N, T — > oo 
with N/T — > q G (0,1], the empirical spectral distribution fiR N weakly converges 
(almost surely) to a probability measure whose density p(x) is defined by: 



1 J (7+ — x)[x — 7_) 
P ^ = 2~^q x 1 b-.7+] dx ( L3 ) 

where j± = 1 + q ± l^fq. 

Independently of the aforementionned work on random matrix theory, much work 
has been devoted to studying the statistics of financial stocks. It turns out that most 
financial assets (stocks, indices, etc..) possess universal features, called stylized facts. 
We refer to the review [5] for a discussion on this topic. On the other hand, many 
models have been proposed in the litterature that take into account these stylized 
facts. Among them, there has been growing interest in the Multifractal Random 
Walk (MRW) model introduced in [1]. The MRW is simply defined as: 

X t = B(M[0,t]) (1.4) 

where B is a standard Brownian motion and M is an independent multifractal 
random measure (MRM for short), see section 2.1 for a reminder of the construc- 
tion/definition as well as (standard) notations used throughout the paper. 

We thus aim at studying the large sample covariance matrices where the un- 
derlyings evolve independenly as Brownian motions with a time change, which can 
be thought of as a volatility process with memory (i.e. the volatility process is cor- 
related in time). The main example of such processes we are interested in is the 
multifractal random walk but we will also consider other examples. More precisely, 
the matrix can be defined, for 1 ^ i ^ N, 1 ^ j ^ T, by: 

x »&) = B m H o4) ~ B mo^ = riU) (1-5) 

where the Bi are i.i.d. Brownian motions and the Mj are i.i.d. multifractal random 
measures independent of the B; L . 

The purpose of this work is to characterize, in the limit of large matrices (N, T — > 
00 with N/T — > q G (0, 1]), the limit of the empirical spectral measure hr n = 
1/-/V"5^i = i $\i where (Ai, . . . , Xn) are the eigenvalues of the empirical covariance ma- 
trix Rn := XjyXjy 

The next sections are organized as follows. In section 2, we remind the definition 
of MRW and introduce the main notations of the paper. In section 3, we state our 
main theorems which are characterizations of the limiting spectral measure of Rn 
through its Stieltjes transform for different type of underlying process X. Since these 
equations are tedious to invert, we leave the analysis of the underlying probability 
measure to a forthcoming study (we are working on this subject). The proofs appear 
in section 4 with some auxiliary lemmas proved in the appendix. The strategy of 
our proofs is classical among the random matrix litterature (the so-called resolvent 
method) as it relies on the Schur recursion formula for the Stieltjes tranform; in 
particular, we follow the approach of [3]. The main difficulty lies in handling the 
Stieltjes transforms in a multifractal setting. 



3 



2. Background, notations and main results 



2.1 Reminder of the construction of MRM 

To fix precisely the notations that we will use throughout the paper, we quickly 
remind the main steps of the construction of Multifractal Random Measures (MRM). 
The description is necessarily concise and the reader is referred to [1] for further 
details. In particular, we use the same notations as in [1] to facilitate the reading. 
We consider the characteristic function of an infinitely divisible random variable Z, 
which can be written as E[e ipZ ] = where (Levy-Khintchine's formula): 

<p(p) = imp - --f 2 p 2 + ! (e ipx - 1) u(dx) (2.1) 
2 Jr* 

and v(dx) is a so-called Levy measure (ie satisfying j Rt min(l, x 2 ) u(dx) < +oo) 
together with the following additional assumption: 



/ \x\ v{dx) < +oo, (2.2) 



so that its characteristic function perfectly makes sense as written in (2.1). We also 
introduce the Laplace exponent ip of Z by ip{p) = <p(—ip) for each p such that both 
terms of the equality make sense, and we assume that the following renormalization 
condition holds: ^(l) = 0. 

We further consider the half-space S = {(t,y);t G R, y G ffi+}, with which we 
associate the measure (on the Borel cr-algebra B(S)): 

6{dt,dy) = y~ 2 dtdy. (2.3) 

Then we consider an independently scattered infinitely divisible random measure \x 
associated to (<p, 6) and distributed on S. 

Then we define a process oj e for e > by the following. Given a positive parameter 
r, let us define the function / : M + — > R by: 



f(r) 



r, if r ^ t 
t if r > t 



The cone-like subset A e (t) of S is defined by: 

A e (t) = {(s,y)ES;y> e, -f(y)/2 ^s-t^ f(y)/2}. (2.4) 

We then define the stationary process (cu e (t))j e R by: 

u e (t)= fi(A e (t)). (2.5) 

The Radon measure M is then defined as the almost sure limit (in the sense of 
weak convergence of Radon measures) by: 



M(A) = lim MJA) = lim / e^ {r) dr 
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for any Lebesgue measurable subset Acl. The convergence is ensured by the fact 
that the family (M e (A)) t>0 is a right-continuous positive martingale. The structure 
exponent of M is defined by: 

Vp ^ o, C(p) = p - 4>(p) 

for all p such that the right-hand side makes sense. The measure M is different from 
if and only if there exists e > such that ((1 + e) > 1, (or equivalently < 1). 

In that case, we have: 

Theorem 2.1. The measure M is stationary and satisfies the exact stochastic 
scale invariance property: for any A G]0, 1], 

(M(XA)) AcB{0>T) l = (Xe^M(A)) AcB{0}T) , 

where Q\ is an infinitely divisible random variable, independent of (M(A))acb(o,t), 
the law of which is characterized by: 

E[e ipnx ] = X~ vip) . 



2.2 Notations 

Let N and T := T(N) be two integers, the aim of this paper is to compute the 
empirical spectral measure of the matrix Rn '■= X^Xn as N — > oo, where Xn is a 
N xT real matrix the entries of which are given by (1.5). Recall that the number iV 
of sampled processes is supposed to be comparable with the sample size T := T(N), 
and more precisely, we will suppose in the following that there exists a parameter 
q e]0, 1] such that: 

N . . 

hm — = q. (2.6) 

iV-»oo I 

We further set Rn '■= ^nXn, and if M is a symmetric real matrix, we will denote 
by ix m the empirical spectral measure of M. 
Define the (T + N) x (T + N) matrix B N by: 

R f X N 

Bn ~{x n 

We also define for z G C \ R, 

A N (z) = (zI T+N - B N ) = 

Notice that 



zI T -X N 
-Xn zIn 



B 



2 _ I fhN 

N 



R N 
R 



N 
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and that the eigenvalues of Rn are those of Rn augmented with T — N zero eigen- 
values. We thus have: 

N T-N. , . 

where 8 X stands for the Dirac mass at x. Combining this equality with the relation 

f(^)fx B 2 N (dx) = J f(x 2 )fiB N (dx) (2.8) 

true for all bounded continuous functions / on R, we see that it is sufficient to study 
the weak convergence of the spectral measure of B N for the study of the convergence 
of the spectral measure /jlr n . 

We will thus work on the (weak) convergence of the spectral measures (Ib n and 
E [fiB N ] in the following. To that purpose, it is sufficient to prove the convergence of 
the Stieltjes transform of these two measures. Recall that, for a probability measure 
(jl on R, the Stieltjes transform G^ of /i is defined, for all z 6 C \ R, as: 

G^z) = [ —fiidx). (2.9) 
J R z — X 

and one can note that: 

G, B Jz) = _l_Trace(G^(z)), (2.10) 

where we have set: 

G N {z) = {A N {z))- x . (2.11) 

Hence, we have to investigate the convergence of the right-hand side of (2.10). Let 
us introduce the two following complex measures L N Z and L 2 N Z such that, for all 
bounded and measurable function / : [0, 1] — > R: 



L N Z (f) - T7 y_]f ( T7 ) G N (z) k+TMT 



T 

y^/ y G N (z) kk 

k=l ^ ' 

k=i v 7 



Clearly, we have the relation 



1 T N 

-Trace(G JV (z)) = — — L^([0,1]) + — — 4 2 ([0,1]) (2.12) 



N + T ' v v " N + T JV VL ' " N + T 
so that it suffices to establish the convergence of the two complex measures L N Z and 

7-2,2 

N ■ 
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3. Main results 



3.1 Lognormal multifractal random walk 

We first present our results when the process X(t) is a lognormal multifractal random 
walk, i.e. X(t) = B(M[0; t]) where M is the MRM whose characteristic and structure 
exponent (see section 2.1) are respectively given by: 

( \ V 7 2 2 

C(g) = (l + y)9-yg 2 - 

We will make the assumption that the intermitency parameter r y 2 is small enough 
so as to overcome in our proofs the strong correlations of the model. 

Assumption 3.1. More precisely, let us suppose that: 

7 2 < ~ (3.1) 

Though we conjecture that our results hold as soon as the measure M is non 
degenerated, i.e. 7 2 < 2 (see [1]), Assumption 3.1 is largely sufficient to cover most 
practical applications. For instance, in financial applications or in the field of tur- 
bulence, 7 2 is found empirically around 2.10 -2 . 

We can now state our main result about the convergence of the empirical spectral 
measures and mean empirical spectral measures of the matrices and Rm'- 

Theorem 3.2. i) There exists a probability measure v onM such that the two mean 
spectral measures K\p BN ] and E[[i Rn ] converge weakly respectively towards the two 
probability measures j^v + o~o and v o (x 2 ) -1 as N goes to oo, where v o (x 2 ) -1 
is the pushforward of the measure v by the mapping x i— )■ x 2 . 

ii) The two spectral measures (jlb n and [ir n converge weakly in probability re- 
spectively to the two probability measures + j^5 and v o (x 2 ) -1 as N goes 
to oo. More precisely, for any bounded and continuous function f , J f(x)fj,R N (dx) 
converges in probability to J f(x)v o (x 2 )~ 1 (dx) . 

Hi) Let Nk be an increasing sequence of integers such that Yl'kLi^k 1 < +oo, 
then the two sequences f^B Nk and Hn Nk converge weakly almost surely to the two 
probability measures j^v + j^S and v o (x 2 ) -1 as k goes to oo. 

Theorem 3.2 is implied by (2.10), (2.12) and by Theorem 3.3: 

Theorem 3.3. %) The measures E[L^ Z ] and E[L^ 2 ] converge weakly towards two 
complex measures. More precisely, there exist a unique [i 2 G C and a unique bounded 
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measurable function K z (x) over [0, 1] such that, for all bounded and continuous func- 
tion f on [0, 1], we have respectively: 



E[L)f{f)] ^n^oo [ K z (x)f{x)dx, 
Jo 

E [!#(/)] ^N^fil [ f(x)dx. 

Jo 



ii) In addition, we have the following relation between /i^gC and K z (x): 

K z (x)dx = qfi 2 z + ^—^ (3.2) 



Hi) Furthermore, there exists a unique probability measure v onl whose Stieltjes 
transform is [i 2 , meaning that for all z £ C \ R ; 

"5 = / r? < 33 > 

It is important to state a characterization of the probability measure v : it is done 
by means of its Stieltjes transform /j, 2 : 

Theorem 3.4. The constant fi z and the bounded function K z (x) are uniquely de- 
termined for all z £ C \ R ; &?/ the following system of equations: 



£ = E 



K z {x) = \ z-qE 



l 

K z {t)M{dt) 



1 ' K z (t)M(dt) 



\t — x\J + 

where M is the MRM with structure exponent ((q) = (1 + r y 2 /2)q — q 2r f 2 /2. 



(3.4) 
(3.5) 



Let us notice that one can give a precise meaning to (3.5) for all j 2 £ [0, 2[. 
Indeed, we can define for all x £ [0, 1] and all continuous function /, the following 
almost sure limit SIS db definition: 

One can also check with this definition that we have: 



i , _ . j2 rl 



f(t)M(dt) = lim / e «w(««(t)^«(*))y( t ) e ««W^ 



\\t — X\J+ ^°J0 

Conjecture 3.5. With this extended definition, we conjecture that theorem 3.4 holds 
in the lognormal multifractal case for all r y 2 £ [0, 2[ and thus that the limiting equa- 
tions can be obtained by the ones of theorem 3.9 (see below) with 2W = u e as e — > 0. 
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3.2 General multifractal random walk 



We now look at the more general case when the change of time is a measure M 
for which the function (f(q) is given by (2.1) and the structure exponent by £(q) = 
q — ip(q) with ip(q) = <p(—iq). 

We still have to make an assumption to avoid the issue of strong correlations. In 
this more general setting, Assumption (3.1) becomes: 

Assumption 3.6. Assume that the structure exponent of the MRM satisfies the 
condition: 

C(2)>5-4C'(1). (3.7) 
and that there exists 5 > such that: 



C(2 + 5) > 1. 



(3. 



As in the previous section, we conjecture that our results hold as soon as the 
measure M is non degenerated, i.e. (see [1]) ((1 + e) > 1 for some e > 0. 

Theorems 3.2 and 3.3 remain unchanged for this more general context. Theorem 
3.4 becomes: 

Theorem 3.7. The constant [i 2 z and the bounded function K z (x) are uniquely de- 
termined for all z £ C \ R ; by the following system of equations: 



E 



: - / K z (t)M(dt) 
o 



K z {x) = [z-qE 



z — 



t 



X 



K z {t)Q{dt) 



(3.9) 
(3.10) 



with k = '0(2) and where M is the MRM whose characteristic and structure expo- 
nent are respectively l p{q),C{ ( l) an d where the random Radon measure Q is defined, 
conditionally on M , as the almost sure weak limit as e goes to of the family of 
random measures Q e (dt) := e Ue ^'M(dt) where, for each e > ; the random process 
ZJ e is independent of M and defined as cu e (t) = ~p(A t (t)) where ~p is the independently 
scattered log infinitely divisible random measure associated to (<p, 6(- fl Aq(x))) with: 



<p{p) = ip^ 2 -k)+ (e ipx - l)(e x - l)v(dx). 



(3.11) 



3.3 Lognormal random walk 

Let us mention that one can easily adapt the methods used to prove the above 
theorems in the simpler case (lognormal case) where X(t) is defined, for all t e [0; 1], 
by: 

X(t) = B ^ e 2W ^ds^j , (3.12) 
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where (W(s)) a6 [o ; i] is a stationary gaussian process with expectation m and station- 
ary covariance kernel k. The normalisation will be chosen such that: m = —k(0). 
In this context, the entries of are given, for 1 ^ i ^ N, l^j^T by: 

1 



X N (ij) 



e Wi(j,) B i ._ 



(3.13) 



where the (Bj)^ are i.i.d standard centered Gaussian random variables and the Wi 
are i.i.d stationary Gaussian processes with expectation m and stationary covariance 
kernel k. Indeed, if one makes the following extra assumption: 

Assumption 3.8. Assume that for some constants C > and (3 > 0, the covariance 
kernel k satisfies: 

VxgM, \k(x) - fc(0)| < c \ x f- 
With the same notations as in the previous section, we can now state the following 
theorem under assumption 3.8: 

Theorem 3.9. The system of equations for [i 2 z and K z (x) becomes: 



K z {x) 



E 



K z (t)e 2W(t Ut 



z-qE 



K z (t)e ik ^e 2W ^dt 



(3.14) 



(3.15) 



where (W(t)) t £[o-i] is a stationary gaussian process with expectation m and stationary 
covariance kernel k. 



4. Proofs of the main results 

In this section, we give the proofs of theorems 3.2, 3.3 and 3.4. The proof of Theorem 
3.7 is very similar and we will not explain it in every detail, except for the final part 
where we establish the second equation of the system in Theorem 3.7 verified by 
K z . We will give the details for this part of the proof in the appendix. The proof 
of theorem 3.9 is an easy adaptation of our proofs for theorems 3.2, 3.3 and 3.4; it 
is left to the reader. Furthermore, the proofs are very similar when q = 1 or when 
q < 1. For the sake of clarity, we assume T = N and hence q = 1 in the proofs that 
follow. 

Hence, in the following, we will suppose (unless otherwise stated) that: 

7 2 2 7 2 
<p(q) = + Q y> 

C(g) = (l + y)g + 9 2 y, 

2 1 

7 < 
' 3' 
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and M will be the MRM whose structure exponent is ( (see section 2.1 for a re- 
minder) . 

Our approach to show the convergence of E[L^] and E[L^ 2 ] consists in proving 
tightness and characterizing uniquely the possible limit points. The classical Schur 
complement formula is our basic linear algebraic tool to study E[L^ Z ] and E[L^ Z ] 
recursively on the dimension N, as is usual when the resolvent method is used. The 
original part of our proof is that we apply the Schur complement formula two times 
in a row to find the second equation of the system in theorem 3.4 involving the limit 
point K z (x) of the measure E[L^]. We will also show that the limit points of the 
two complex measures E[L^ 2 ] and E[L^ Z ] satisfy a fixed point system (written in 
theorem 3.4). 

We begin by showing tightness. 

4.1 Tightness of the complex measures E[LJ],E[LJ] and 
limit points 

Lemma 4.1. The two families of complex measures (E[L^])jveN, i = 1,2 are tight 
and bounded in total variation. 

Proof. Let us present the proof for (E[L^ z ])7veN; the other proof is similar. 
One has, for each iV G N: 

i N i 
I E[2#] | [0, 1] = -^ I nG N (z) kk ] | < —— (4.1) 

k=l \ \ ) \ 

and so the family of complex measures (E[L]^ z ])jv e N is bounded in total variation. 
It is obviously tight since the support of all the complex measures in the family is 
included in [0, 1], which is a compact set. □ 
Using Prokhorov's theorem, we know that those two families of complex measures 
are sequentially compact in the space of complex Borel measure on [0, 1] equipped 
with the topology of weak convergence. In particular, there exists a subsequence 
such that, for all bounded continuous function /, one has, when iV goes to +oo 
along this subsequence: 

E [!#(/)] / f(x)nl(dx). (4.2) 
J o 

Lemma 4.2. The complex measure y} z (dx) has Lebesgue density; more precisely, 
there exists a bounded measurable function K z (x) such that: 

nl(dx) = K z (x)dx. (4.3) 
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Proof. One has: 



1 N 

E [/#(/)] | < -J2\f( k / N )\ E l G ^)kk] (4.4) 



fc=i 

A? 



Letting iV — » +oo along a subsequence, one obtains: 



f(x)j/ z (dx) 



1 



< T^y I \fW\ dx - ( 4 - 6 ) 



This proves the lemma. □ 
Thus, there exists a subsequence such that, as N tends to +oo along this subse- 
quence: 

E [L]f(f)] / f(x)K z (x)dx. (4.7) 
Jo 

Lemma 4.3. There exists a subsequence and a constant ji 2 z G C such that, as N 
goes to +oo along this subsequence: 

E[L% z (f)] ! f(x)dx. (4.8) 

Jo 

Proof. It is easy to see that the G]y(z) kk ,k = N+l, . . . , N are identically distributed. 
In particular, these variables have the same mean (j%(N). One has, for all N: 

< 7c*7I\T ■ ( 4 -9) 



So there exists a subsequence and a complex number p? z such that, as N goes to +oo 
along this subsequence, fJ%(N) — > fi 2 z . One thus obtains, as N goes to +oo along this 
subsequence: 

E[L% z (f)} -+fif f(x)dx. (4.10) 
Jo 

□ 

Following the classical method as in [3], [2], [7], we will show in the following 
that the limit point (i 2 and K z (x) are defined uniquely and do not depend on the 
subsequence. We will first recall some preliminary results on resolvents. 



4.2 Preliminary results on resolvents 

We first recall the following standard and general result; the next lemmas of this 
section are also standard but are applied to our particular case. 



12 



Lemma 4.4. Let A be a symmetric real valued matrix of size N. For zeC\l, let 
us denote by G(z) the matrix 

G{z) = (z-A)- 1 . (4.11) 
For z£C\K and k G {1, . . . , N}, we have 

S(z)3(G(z) fcfc ) < and |G(z) fcfc | ^ (4.12) 

In particular, if F C {1, . . . , iV} z's a finite set and (ai)i e p a finite sequence of positive 
number, then: 



and we also have: 



^ 1. (4.13) 



p - Eiejr aiG(z)u | |9f(z) 

Proof. Write A = U t DU where D is a diagonal matrix with diagonal real entries 
{h)x^i^N- Then 

- 1 

G{z) kk = Y^\U ki \ 2 

Z — A; 



i=l 



Since = (S( 2 )-A0*+3( 2 ) a and ^(A) = pR^W the elation (4.12) 

follows. It is then straightforward to derive (4.13) from (4.12). □ 
For i = 1, . . . , N, let Xjj = (X^(kl)) k ^ be the matrix obtained from Xn by 
taking off the z-th column and row. Define, also for % = 1, . . . , 2N the (2N — 1) x 
(2N — 1) matrix A$(z) obtained from A^(z) by taking off the i-th column and row. 
In particular, for i — 1, . . . , N, 



A^ +i \z) - ' zIn 



Fori = l,...,2iV, set: 



X$ zIn- 



G { S{z) = {A^{z))-\ (1.10) 
^ N uenuie me mtiirix X^ 



Let now X$ denote the matrix Xm with the i-th column and row set to and 



Apf(z) denote the matrix Ajy(z) with the i-th column and row set to excepted the 
diagonal term. Again we have, for % = 1, . . . , N: 

N {Z) -\-X$ zI N 
For i = 1, .. .,2N, set: 

&$(z) = {A®(z))-\ (4.16) 
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In the paper, we will also use the terms A^' l \z), G^' l \z), A^' l \z), A^' l \z). The 
double superscript just means that you make the operations described above to the 
rows and columns % and k. 



Lemma 4.5. For all k € {1, . . . , N} and all t ^ N + k, one has: 

1 



E 



Gn(z)u — G^ +k \z) tt 



N\%(z^ 2 



Proof. Multiply the identity: 

A$ +k \z) - A N (z) = A^ +k) (0) - A N (0) 
to the left by Gn(z) and to the right by G^ +k \z) to obtain 

G N (z) - G^ +k \z) = G N (z)(A^ +k \0) - A N (0))G^ +k \z 
Then one has: 

G N {z) tt - G^ +k \z) tt = (G N (z)(A^ +k) (0) - A N (0))G% +k \z) 

N 

= G^ +k \z) N+kfi s ^2G N (z) ti r k (i) 



N 



+ G N (z) t ,N+k ^2r k {j)G { x +k \z) jt 

N 

= G N (z) ttN+k ^2r k (j)G 



where we have noticed that, for all t ^ N + k, G^ +k \z) N+k,t = 0. 
Therefore, we find that: 



E 



G N {z) tt - G$ +k \z) tt ^ E [\G N (z) t , N+k \Y'^ 



21 1/2 



N 



(4.17) 

(4.18) 

(4.19) 

(4.20) 
(4.21) 

(4.22) 
(4.23) 



J2r k U)G ( N N+k \z) Jt 



1/2 



.24) 



by Cauchy-Schwartz's inequality. Using then the independence of r k (j) and G^ +k \z), 
we get: 



E 



G N (z) a - G^ +k \z) tt \] ^ E [\G N (z) t , N+k \ 2 ] 1/2 ^ h(l) 2 ] 1/2 E 



N 



V 



y/N\S(z)\ 2 ' 



The proof is complete. 



□ 
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Lemma 4.6. There exists a constant C > such that, for all k e {1, . . . , N} and 
all t 7^ k: 



E 



G N (z) tt -G%\z) tt 



C 



1 



3(*)I 2 JV^' 



(4.25) 



Proof. Again, we start from the relation: 

G N (z) - G%\z) = G N {z)(A%\0) - A N (0))&£\z). 

Thus we have 

G N (z) tt - &*\z)u = (g n (z)(A%\0) - A N (0))G%\zj) 



N 



= G%\z) k , t G N (z) u n(k) 

i=N+l 
N+l 

+ G N (z) ttk J £r j (k)G ( £\z) jt 



3=1 
N+l 

= G N (z) tjk J2rj(k)G 

3=1 

where we have noticed that, for all t ^ k, G^\z) k ,t = 0. 
Therefore, we find that: 



N \ Z >3t 



(4.26) 
(4.27) 

(4.28) 

(4.29) 



E 



G N (z) tt - &*\z) u < E [\G N (z) ttk \ 2 ] 1,z E 



21 1/2 



N 



-| 1/2 



(4.30) 



J2r,{k)GP(zh 

3=1 

by Cauchy-Schwartz's inequality. We want to expand the square in the above expres- 
sion. To that purpose, we first observe that, conditionnally to the M l , the variables 
(rj(k))j are independent from G$(z) and centered. Hence we have for j ^ j', 



E 



rj (k)r f (k)G%\z) Jt G%>(z) ft 



(*)/ 



0. 



Thus we get: 



E 



G N {z) tt -G^\z) tt ^E[\G N (z) tik \ 2 ] 1/2 [ J2 E 



'N+l 



^E[|G w (^, fe | 2 ] 1/2 (E E b( A; ) 4 ] 1/2E 

J=i 



J=i 

'N+l 



r 3 {kf 



G^\z) jt 



1/2 



1/2^ 



1/2 



Ejr^A;) 4 ] 1 / 4 
E[n(A:) 4 ] 1/4 



''N+l 



1/2^ 



1/2 



(JV+1 
3=1 



\G%\z)j t 



1/4 
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Now we use the scaling properties of the MRM to obtain, for some positive constant 
C, 

E[ rj (k) 4 } = 3E[M(0, -^) 2 ] ^ CN-« 2 \ 
Furthermore, by using Lemma A.l which assures that, almost surely: 



N+l 

El<* 



2 1 



\%(z)\ 2 



(4.31] 



and the fact that: 



we finally obtain 



E 



JV+l 



G N (z 



f N+l 



\ G N » - './ 



( z )jt 



N\Z tt 



(4.32) 



|3(z)| 2V iV 



□ 



(4.33) 



r N \*Jtt 

It just remains to check that £(2) = 2 — r y 2 . 
Lemma 4.7. For each k £ {1, . . . ,2iV} ; ift^k, then 

G£ (z) tt = G^ N \z)u, 
and ift — k, then G^\z) kyk = z~ x . 

Proof. It is straigthforward to see that the two matrices G£\z) and G^\z) have the 
same eigenvalues except that G^ (z) has one more zero eigenvalue. In addition, the 
eigenvectors look also very similar since you can obtain 2N eigenvectors of G^\z) 
by adding a zero entry to the eigenvectors of G^ (z) (between the entries k — 1 and 
k). The last eigenvector of G^ \z) is the vector of ~R N for which all entries are zero 
except the entry number k. 

Now observe that with G%\z) = Udiag(z-X)U* and G%\z) = Kdiag(z - X)V*, 



2N 



Gm (z)u — ^2 



1 



i=i 

N 



gP(z)u=y: 



Uti 



Vti 



i=l 



z - Xi 

■ 1 

z - Xi 



(4.34) 
(4.35) 



The result follows since, for t ^ k, 

2N-1 



E 

i=i 



z - Xi 



2N 

E 

i=l 



Vti 



z - Xi 



and, for t = k,G^\z)k, 



(4.36) 
□ 
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Lemma 4.8. For all z G C and Lebesgue almost every point x G [0, 1], we have 

S(z)$r(K z (x)) ^ (4.37) 

and 

|9(tf,(:n))| < -1- (4.38) 

Proof. This is a straightforward consequence of Lemma 4.4. Indeed, we have for all 
positive continuous function / on [0, 1] and N G N: 



We pass to the limit as N goes to oo along some suitable subsequence and obtain: 



%{z)%(J^ f(x)K z (x)dx^ ^0. 



The result follows. □ 

4.3 Concentration inequalities 

This lemma is adpated to our case from Lemma 5.4 in [3]. 

Lemma 4.9. Let f : [0, 1] — > R be a bounded measurable function. For each i G 
{1,2} ; we have the following concentration results: 

E [I !#(/) - E[I#(/)] | 2 ] < ^^Jy • (4-39) 

Proof. Define two functions F^ and such that: 

Fi ( , < , < ™ • 1 < < « N ) = h t f (I) ("») 

/ \ 1 N+l / U \ 

^{{ x S n ) ltltK+1 ' l<i<N ) =Fg / (ivTT) GKWw+K (441) 

We will prove the Lemma for L]^ z ; the proof for L 2 ^ is a straightforward adap- 
tation. 

Let, for k G {1, . . . , JV + 1}, 



^=^((^r) 1<i<iV . ( 4 - 42 ) 
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If P denotes the law of the vector [x[^ j 
E[|ift-E[ift] | 2 ] 

N 

= j>0 E K i - E[i^ i n i 2 ] 

i=0 

= E/ ' / (^v( x l» x 2,---,Z»+l,yi+2,---,^ I 2 
i=0 ^ ^ 

ciP 0l+1 (x) 

, x 2 , . . . , Xj, x i+2 , xat + i) - F N (xi,x 2 , ...,Xi,y, x i+2 , x N+1 )) dP(y) \ 

i=0 

dP m+1 (x) 

N 

^ SU P II V^Fjv || 2 / || 2 dP® 2 (x,y). 

i=0 R(iv+i) 2 J 

The quantity V Xj+1 Fjy refers to the gradient of in the direction of the vector 

If we consider a couple of processes (.B^M 1 ) independent from (P^M 1 ) with 
the same law, it is easy to see that: 

r N 

/ || x-y || 2 dP®dP(x,y) = E E H 1( o-) -^(o,^) -^(o-) +^(o,^)) 2 ]- 

J 3=1 

N 

= 2 - 2j2 E [( Bl M i(o4) ~ S lfi(o,^))(^l?i(o ±)^h(o,i^y 

3=1 

= 2. 

In our case, we have, for i e {1, . . . , TV + 1}, j e {1, . . . , iV}: 

= GN{z)k,jGN{z)N+i,k + GN{z)k,N+iGN{z)j,k (4.43) 

Thus, 

V, l+1 Pv = ^ E / ( liF ) V, l+1 G^(^) fcfc (4.44) 
k=i ^ ' 

It is now plain to compute: 



V Xi+1 Fjv 



1 * 

2 - -£| (G JV (z) J D 1 (/)G' J v(z)) Ar+mi .+ (G' J v(z) J D 1 (/)G'^(z)) : 



i=i 
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where D 1 (f) is the (2iV)-dimensional diagonal matrix of entries: 

D 1 (f)kk = f (lyj -k 1 < k < N }- 



One thus has: 



N 



1 ||2 



3=1 
2N 



3=1 



4 


1 f II 2 

\ J Woo 


AT2 


%Z | 4 



where, in the last line, we used lemma A.l and the fact that the matrix G^{z)D x {f)G^{z) 
has a spectral radius smaller than || / ||oo / | ^ z | 2 - 



Finally, 



E[|F* -E[F£] I 2 ] < 



8 1 




2 

oo 


N 


Qz 


1 



(4.45) 

□ 



We also prove the following lemma: 
Lemma 4.10. For all a > 1 such that C(2ct) > 1, we have 

N 

5>(*) a (Gr fc) w«-E[Gr fe) (-)« 



E 



t=i 



CQniV) 5 



(4.46) 



/or some positive constant C independent from N,z,k. 

Proof. Notice that (rk(t)) t and (z) are independent. Hence, by conditioning 

with respect to the process (rk(t)) t , we can argue along the same lines as in the 
previous lemma with rk(t) instead of ■hf(4f) and we get the formula: 



E 



N 

J2r k (tf (&» + %)u-E[&» + %)u] 
t=i 



\S(z)\± 



E[supr fc (t) 4 ]. 
t 



We conclude with Proposition B.l in the appendix . □ 
In the following, we fix a > 1 such that C(2at) > 1 (because of the expression of 
( and the inequality ■y 2 < 1/3, it is clear that such a number exists). 
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4.4 The system verified by the limit point yr z and K z (x): first 
equation 

From the Schur complement formula (see e.g. Lemma 4.2 in [3] for a reminder), one 
has for k G {1,...,N}: 



GN{Z)N+k,N+k — 

Using Lemma A. 2, one can write: 



N 



Z — 



St 



,t=l 



-1 



(4.47) 



G N (z 



Ny^)N+k,N+k 



N 



z-Y,r k {tfG^ +k \z) tt + e^ k {z) 
t=i 



(4.48) 



where e] Vfe (^) is a complex valued random variable for which there exists C > such 
that for all N G N and 1 < k < N, 



C 



(4.49) 



By using Lemma 4.7, we can write: 

N 



Gn(z 



N\ z )N+k,N+k 



z)tt + e Nh (z) 



t=l 



(4.50) 



Lemma 4.10 applied to a > 1 such that ({2a) > 1 yields: 



E 



;(N+k) , 



J2rk(t) 2 (G^ +k \z) tt -E[G^ N y , m 
Thus, one can write: 



C(lniV) ; 



(4.51) 



Gm(Z 



N{ z )N+k,N+k 



N 



$>,(t) 2 E G% +k \z) tt 



t=l 



(4.52) 



where £% k {z) * s a complex valued random variable such that for all iV G N and 

(4.53) 
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In addition, using Lemma 4.5, we can show: 



E 



N 



t=i 

N 



< ^E[r fc (t) 2 ]E 



G^ +k \z) tt -G N (z) tt 
G$ +k \z)u-G N (z)tt 



t=i 



1 



\&(z)\WN 



(4.54) 

(4.55) 
(4.56) 



It follows: 

GN{z)N+k,N+h 



N 



^r fe (t) 2 E [G N {z)u\ + e x Ntk (z) + e\ k {z) + 



t=i 



(4.57) 

where € 3 Nk (z) is a complex valued random variable such that for all iV 6 N and 
1 ^ k < Ar+ l, 

E[l&,fc(*)l] < ^ i ^ - (4-58) 



|3(z)|VW 
hei 

Lemma 4.11. TTie following inequality holds: 



Let us denote by the interval [Kf, ^]. Then we have: 



E 



A' 



^(r^f-M^jE^)* 



t=i 



/or some positive constant C . 

Proof. We expand the square and, because r k (t) and rk(t') are independent for t ^ t! 
conditionally to M k , we have: 

N 



t=l 



E\\^[r k (t) 2 - M\P N ))E[G N (z) a ] 

N 

£ E[(r*(t) a - M fc (4)) (r fc (t') 2 - M fc (/^) )e [G N (z) tt ] E [G,( Z ) W ] 
t,f=i 

AT 2 

^E[(r fc (t) 2 -M fc (4)) ]e[G,v^)«] 2 
t=i 

2V 

2^E[(M fc (4)) 2 ]E[G iV (^ t ] 2 



£ 2C 



N 



MW\$r(z)\' 
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We can thus write 

GN(z)N+k,N+k 



N 



£m^)E[GW(*)«] 



+ ejf Jz) + e% h (z) + h (z) + e% h {z 



-N,k 



-1 



□ 

(4.59) 
(4.60) 



where k {z) is a complex valued random variable such that for all N G N and 

Set e N;k (z) = e l Nk {z) + e 2 Nk (z) + e 3 N Jz) + e% Jz) and rewrite: 



G N {z 



N\Z)N+k,N+k 



N 



J2M k (I t N )E[G N (z) tt ] + e N , k (z) 



t=i 



-i 



(4.62) 



We now need to introduce the truncated Radon measure M k (dx) with Lebesgue 
density e Ue ^ which converges almost surely as e goes to 0, in the sense of weak 
convergence in the space of Radon measure, to the measure M k (see section 2.1). 

Lemma 4.12. For e > 0, the following uniform bound holds: 



N 



N 



mvE[\ y £ / M k (I t N )E[G N (z)tt] - ]T M £ fc (4)E [G N (z) tt ] | 2 ] < 

N + i + 1 



t=i 



WzW' 



Proof. We expand the square. Note that the covariance function p e of the process 
u e increases as e decreases to and uniformly converges as e — > towards ln + o 
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t=l 



over the complement of any ball centered at 0. Thus we have: 

N N 

supE[| Y,M k (I l N )E [G N {z) tt ] -J2 M eK)® [G N (z) tt ] | 2 ] 

t=i 

AT 

snp ^ E[(M fc (4) -M*(4))(M fc (iS) - M £ fc (4))]E E 
N t,t'=i 

TV 

snp £ E [(M fc (4) - M e fe (4))(M fe (4) - M e fc (4))j E [G N (z) tt ] E [G, 



'N\Z)t't' 



t,t'=l 
N 



= sup V E M fc (4)M fc (4) -E M e fc (4)M e fc (4) J E [Gw(z) tt ] E [G N (z) t ' t >] 



t,*'=l 
JV 



= sup V E E [GW(zV,] / / (e Wbl+ ^ - e^^) rfrrfw 

w)| 2 Jo Jo 



where, in the fourth line, we used the fact that, if J-" e is the sigma field generated by 
the random variables fJ,(A), A E B({(t,y) : y ^ e}), then E[M k (A)\T £ ] = M*{A) for 
all borelian set A. A straightforward computation leads to the relation 



ln^ + l-| if|t|<e 
Pe(t) = { lug ' ife^|t|<r 

ifr<|t| 

By using the expression of p e , it is then plain to obtain the desired bound. 
We can thus write 



N 



N+k,N+k - 

where 
and also: 

E [Gjv(z)jV+fc,7V+fc] = E 



M e fe (4)E [G N (z) tt ] + e N , k (z) + 8(e, N, 



t=i 



supE[|<5(e,A^,z)| 2 ] ->■ as e ->■ 0, 

N 



(4.63) 



□ 



(4.64) 



(4.65) 



N 



- M*(I l N )E [G N (z) tt ] + e N , k (z) + 5(e, N, . 



t=\ 



(4.66) 

The next step is to study the convergence of the above quantity. Hence we prove 
(see the proof in the appendix): 
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Lemma 4.13. The random variable YltLi M^(Ij^)K [Gn{z)u] converges in probabil- 
ity as N +00 towards K z (x)M k (dx). 

We fix e > 0. For that e, the family of random variables (S(e, N, z))n is bounded 
in L 2 so that it is tight. Even if it means extracting again a subsequence we assume 
that the couple ($^t=i M^I^W, [Gn(z) u ] , 5(e, N, z))n converges in law towards the 
couple (J^ K z (x)M k (dx),Y e ). We remind the reader of (4.47) which implies that 



N 

- M^I^E [G N (z) tt ] + e N>k (z) + 5(e, N, z) 



The quantity [z - E t =i M k (P N )E [G N (z) tt ] + e N>k (z) + S(e, N, z) 
uniformly with respect to N, e and converges in law towards 

-1 



is therefore bounded 



z- j\ z {x)M k t (dx) +Y)j 



We deduce that the expectation of the former quantity converges as e — > towards 
the expectation of the latter quantity. From (4.66), we deduce that 



l'~z 



E 



z- / K z {x)M k e {dx) + Y t 



-r 



(4.67) 



Clearly, standard arguments prove that J Q K z (x)M k (dx) converges almost surely 
towards K z (x)M k (dx) as e — >■ [K z is deterministic (see lemma 4.9), measurable 
and bounded) and, because of (4.65), Y e converges almost surely towards as e — > 0. 

Again, because the quantity (z — K z {x)M k {dx) + Y^j is bounded uniformly 
with respect to e, we deduce that: 



i4 



E 



K z {x)M k {dx\ 



(4.68) 



4.5 Second equation 

Now we turn our attention to the terms Gjv(z)fcfc for k G {1, . . . , N}. Again, by using 
the Schur complement formula, we can write, for k G {1, . . . , N}: 

-1 

(4.69) 

-1 

(4.70) 



N 



GN(z)kk 



J2^(k) rj (k)G^\z) N+i. 



N+j 



N 



^r i (k) 2 G { ^\z) N+l)N+i + 77^ 
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where, using Lemma A. 3, r)j^k(z) is a complex valued random variable for which 
there exists c> such that for all iV e N and 1 < k ^ N, E[|r^ fc (,z)| 2 ] < c/N. 

With a further use of the Schur complement formula for the term G^\z)N+i,N+i, 
we obtain: 



G N (z) 



kk 



N 



N 



r,(s) ri (t)G^ N+l \z) 



St 



s,t^k 



-T -1 



(4.71) 



where G%'"^ v (z) = A^^'iz)" 1 . Note that G K ^ v {z) is independent of {ri{t)) t=l _ N . 



?% N+i \z) = Ap N+i \z)-\ Note that Gp N+i) < 
Using the same arguments as in the derivation of the first equation (in particular 
Lemmas A. 2, 4.7, 4.10, B.l, 4.6 and 4.5), one can show that: 



G N (z 



N\Z kk 



N 



-E 



tt z ~Eti M i (I t N )E [G N (z) tt ] + 5 NAl (z) 
where (5jv,fc,i(^))i <; i «; n are complex random variable such that 

C 



-l 



(4.72) 



n\5 N .kAz)\\ ^ 



■ l\-~< 2 C(2a)-1\ 



(4.73) 



for some positive constant C that does not depend on i,N and for a > 1 such that 
C(2a) > 1. 

Lemma 4.14. One can write: 



GN{z)kk 



N 

E 



where 7f N k (z) is a random variable that tends to m probability as N goes to oo 
Proof. By using Lemma 4.4, we deduce that: 



+ VnA z ) + v 2 nA z ) 



(4.74) 



N 

E 



r;(/c) s 



r-i(/c) 2 



i 1 ^ - Eti [G w (z)«] + ^(z) 2 - Eti M i (J* r )E 



mz 



l w 

— - ^r i (A;) 2 min(|5 7V ,fe,i(^)|,2). 



(4.75) 



i=l 



We stress that the lemma is proved as soon as we can prove that the left-hand side 
in (4.75) converges in probability to 0. Hence it is enough to prove that 



E 



N 



22 r i( k ) 2 min(|<W,i(20|, 2 ) 



i=l 
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converges to as N tends to oo. By noticing that: 

N N 

<W*) = E nisy^Gp^izU -^Mt/^E^^] , (4.76) 

s,t^k t=l 

it is straightforward to see that the variables ( rj(fc) 2 min(|^7v k i{ z )\i 2) ) are 

V ' ' / 1 < i ^ N+l 

identically distributed. Thus we have 

N 

E[J2ri(k) 2 min{\S NAi (z)\,2)] = NE[ n (k) 2 mm(\5 N , k>1 {z)\, 2)] . 

i=l 

Then for all A > 1 and a > 0, we have 

WE[n(A;) 2 nim(|£ w (*)|,2)] =AHE [r^/c) 2 imn(\5 Ntkt i(z)\,2)t 

{Nri(k) 2 s: A} J 

+ NE [r 1 (k) 2 mm(\5 N:kA (z)\,2)l {Nri{k) 2 >A} ] 
«C AE[\5 NAl {z)\]+2E[Nr l {k) 2 l {Nri{k) 2 >A} ] 

=-^- + 2 J^1e\mHQ, -r +i ] 

By using the scale invariance property of the measure M 1 , we have: 
in such a way that 

AC Ni'{ 1 + a ) 
7VE[r 1 (A;) 2 min(|^ M (z)|,2)] < + 2E [M\0, l) a+1 ] ——. (4.77) 

N~ — ^ 

Since C(2) > 5 - 4C'(1) (this inequality is clear with C{q) = (1 + l 2 /2)q + gV/2 
and is due to our hypotheses of Assumption 3.7 in the more general case), we can 
choose p > such that 

^ r ± >p> l-C(l) = ij'(l). (4.78) 

The mapping a e]0, +oo[i->- pa — ip(l + a) reduces to for a = and, because 
p > ip'(l), is strictly positive for a > small enough. So we choose a < 1 such that 
pa — t/)(1 + a) > and we set A = iV p . We obtain: 



NE[ ri {k)'wm{\5 Ntk M\,2)] < + 2*^^(0, T)^] ^ (1+a) . 



The result follows by letting iV — > oo since min((£(2) — l)/4 — p, — ^(1 + a)) > 
0. □ 
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Lemma 4.15. There exists a constant c > 0, which does not depend on N, such 
that for each N G N: 



E 



N 



YltLi M i (I t N )E [G N (z) tt ] 



— E 



Proof. The proof is straightforward using the fact that for i e {1, . . . , N}, the 
random variables 

(4.79) 



z-ES=i M i (I t N )E[G N (z) tt 
are i.i.d. random variables and Lemma 4.4. 
Therefore we can write 



□ 



G N (z 



N\Z kk 



N 



-E 



ri{kf 



+ VN,k{ Z ) +V 2 N,k( Z ) +VN,k( Z ) 



(4.80) 



with E[( V ^ k (z)) 2 } < 

Now we can take the expectation in (4.80) to obtain 







k=l 


1 N r 
= ^f(k/N)E 

k=i 




1 " r 
= ^J2f(k/N)E 

k=l 





/V 



r;(/c)' 



+ VN,k( Z ) 



-EtiM(/^)E[G w (^]J 



+ ?7iv,fe(^) 



with r] Njk (z) = Vn k( z ) Jr VN( z ) J '' r }N k( z ) ■ Then, by introducing the truncated measure 
M e and by using the Girsanov formula, we can approximate (uniformly in N) this 
last expression by: 



v 



N 



J2f(k/N)E 



z-NE 



M \—- -1 

lvl e [ N , N \ 



k=l 



-ELM^ElG^zh] 



+ S(N,k,z,e) (4.81) 



with swp Nk E[\6~(N, k, z, e)| 2 ] going to when e is going to 0. Along some appropriate 
subsequence, this latter quantity converges as iV — > +oo to: 



Jo 



f(x)E 



z-E 



U e (x) 



£ K z (r) M e (dr) 



-i 



dx 



(4.82) 
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where Y e is such that E[(y e ) 2 ] converges to when e is going to 0. And, we thus 
obtain gathering the above arguments that: 









f f{x)E 






Jo 




z 



uj e {x) 



K z (r) M £ (dr 



dx. 
(4.83) 



It remains to pass to the limit as e — > in that expression. This job is carried out 
with the help of a Girsanov type transform in Appendix C. □ 



4.6 Uniqueness of the solution to the system of equations 

Let X be the space of bounded measurable functions [0, 1] — > C endowed with the 
uniform norm defined for / G X by: 



oo = SUp |/(x)|. 

x6[0,l] 



(4.84) 



Define the operator T : X — > X by setting, for g 6 X and for all x G [0, 1]: 

1 



Tg(x) 



z-qE 



For g,h G X and for all x G [0, 1], we have: 



\Tg(x) -Th(x)\ ^ 



i / T \r 



E 



\t — x\ 



1 / T X7 



\t — X\ 



g(t) -h(t)\M(dt) 
9 - hWoo 



i / T ^r 



|SU)| 4 ./ n V I* — arl 



M(dt) 
dt\ \g — h\ 



Recall that 7 2 < 1/3, and thus it is easy to see that: 



i / T \r 



sup 

cce[o,l] Jo \ 1^ — x \ 



dt < +oo 



And we can deduce that there exists a positive constant C such that: 



sup \Tg(x) — Th(x)\ ^ 
xe[o,i] 



C 



(4.85) 



(4.86) 



(4.87) 
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If z is such that C /\^s{z)\ A < 1, the operator T is contracting and thus has a unique 
fixed point g in the Banach space X. We conclude that, for each z with large 
enough, there exists a unique bounded function K z : [0, 1] — > C such that for all 
are [0,1]: 

K z (x) = f — . (4.88) 



z - qE 



Using the first equation, it is now plain to see that, for z such that C/|S(z)| 4 < 1, 
the constant p? z is uniquely defined by the system of equations (by the first equation, 
it is a function of the function K z , which is uniquely defined for such z). 

Now it remains to show that the limit point /j, 2 is uniquely defined for all z G C\R. 
It will be easy to see using analycity arguments. Indeed, from the Montel theorem, 
every limit point fi 2 is holomorphic on the set C\R since it is the pointwise limit of a 
subsequence of the sequence of holomorphic functions L]^ z ([0, 1]) that are uniformly 
bounded on each compact set of C\R (see Lemma 4.4). Thus, fjL 2 is uniquely defined 
for each z e C\R by analytic extension (we have just seen that fi 2 is uniquely defined 
for a set of z with accumulation points). 

The same argument holds for the unicity of the integral K z (x)dx. Indeed, 
every limit point K z (x)dx is a holomorphic function on C \ R that has some 
prescribed value on the set {z e C \ R : C/|S(2;)| 4 < 1}, which has accumulation 
points. 



4.7 Proof of Theorem 3.2, 3.3 and 3.4 

Let us gather the above arguments to prove the main theorems. 

Proof of theorem 3.4'- it is a direct consequence of sections 4.4, 4.5 and 4.6. 

Proof of theorem 3.3 i): The limit points K z (x)dx and n 2 z dx of the two complex 
measures E[X]\} 2 ] and E[L^ 2 ] are uniquely defined because [i 2 and K z (x) satisfy a 
fixed point system of equations (we have just seen this in theorem 3.4). 

Proof of theorem 3.3 Hi): We need to prove that /j, 2 is the Stieltjes transform 
of a probability measure v. From [6], it suffices to prove that fi 2 is holomorphic 
over C \ R, maps {z e C \ R; 3f(z) < 0} to {z e C \ R; 9f(z) > 0} and that 
lim.y_j.jxj iy^i 2 y — 1 (y e R). Let us check those properties. We have already seen in 
section 4.6 that jx 2 is holomorphic. From Lemma 4.4, (i 2 maps {z e C\R; ^s(z) < 0} 
to {z G C \ R; ^s(z) > 0}. Finally, from Theorem 3.4, we have 

1 

1 -z- 1 jl K z (x)M(dx) 

As |K z (x)| ^ the term J Q K z (x) M{dx)/ z converges pointwise towards 

when z = iy and y — » oo. Furthermore, from Lemma 4.8, we have Qf(z)$t(K z (x) ^ 
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-1 

in such a way that 



- JlK z {x)M{dx) 



^ 1 - Therefore 

' C 1 



z - J 1 K z (x) M(dx 

when z takes on the form z — iy (y 6 R). The dominated convergence theorem 
then implies that lim^oo iy/i 2 y = 1 and we can conclude /i 2 is indeed the Stieltjes 
transform of a (unique) probability measure v. 

Proof of theorem 3.2 i) and 3.3 ii) We observe that, for z G C \ R: 

z! T \ (z 2 I T - t X N X N -z l X N 



Let us rewrite the matrix Gn(z) = An(z)" 1 under the form: 

G ~<-~> - few • (4 - 90) 

where Gi(z), Gi^z), G 2 (z) are respectively of size TxT,NxT,NxN. 
By taking the inverse in the relation (4.89), we obtain: 



I T /z 0\fG l (z) fc? lj2 (z)\ = f(z 2 I T - t X N X N y 1 B 
-X N /z 2 I N /z \G lt2 (z) G 2 (z) J 1 J n/z 2 



(4.91) 



where B = (z 2 I T - X N X N )- U X N / z. 

It can be rewritten, using the fact that — XnG\{z)+zG\^{z) = and — Xn 1 Gi ) 2{z)+ 
zG 2 {z) = In, as: 



G 1 (z)/z 1 G la {z)/z\ f(z 2 I T - t X N X N )- 1 B 
I N /z 2 J 1 J n/z' 



(4.92) 



Therefore, taking the trace we get: 



T 

— G N (z) kk = ^tr(z 2 I T - XnXnY 1 , (4.93) 

fc=l 

and, by using the fact that the eigenvalues of X N X N are those of X N X N aug- 
mented with T — N zeros: 

1 1 T — N 

jr E G ^ z ^ = y tr ^ 2/jv - XxXxY 1 + "jTi- ( 4 - 94 ) 

k=l 

Now, taking expectation and using theorem 3.3, we deduce: 

/ K z (x)dx = qz lim — E rtr(z 2 i"y - X^-X^r 1 ] + — ^ (4.95) 

7 Af-s>oo N z 
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Using the fact that (by (2.7)) the spectrum of contains 2N eigenvalues which 
are the positive and negative square-roots of the spectrum oi R N = X N X N plus 
T — N zero eigenvalues and the fact that l/(z — A) + l/(z + A) = 2zj [z 2 — A 2 ), we 
can see that: 

— ^ J2G N (z) kk = j^tr(z 2 I N - X N X N r + (4.96) 

k=l 

Using the relation 2.12 and theorem 3.3, it is easy to see that: 

i N+T i / r 1 \ 

J™ WT^ E nG N (z) kk ] = -—[q£ + K z (x)dx (4.97) 
yv^+oo N + T ^ 1 + q \ J J 

Taking expectation in 4.96 and using (4.97), we get: 

1 (,..2 . f 1 rs _ ^ 1 



1 + q 



(qui + [ K z (x)dx] = p^- lim i-E [tr^/jv - X^X^ 1 } (4.98) 
V Jo / 1 + Q N ^°° ^ 

(4.99) 



1-ql 



1 + q z' 

From equations (4.95) and (4.98), we get the following relation: 



1 l-o 
K z (x)dx = qfi 2 z + -. (4.100) 



and theorem 3.3 ii). is proved. 
With (4.100), (4.97) becomes: 



1 1 / 1 — \ 

W 5S» FTt E «P»M»I = TT^ ( 2 « + -T 1 ) (4101) 



and, we note that the right hand side of (4.101) is the Stieltjes transform of the 
measure 2q/ (1 + q)v(dx) + (1 — g)/(l + q)5 (dx). Thus, the mean spectral measure 
E[/i Bjv ] converges weakly to the measure 2g/(l + q)v(dx) + (1 — g)/(l + q)8 (dx). 
We have also: 

lim -j-E \tr{z 2 I N - XnXn)- 1 } = ^ (4.102) 

Again using the fact that, for all iGl, 1/(z 2 — x 2 ) = (1/ (z — x) + 1/ (z +x)) / (2z) 
and the fact that v(dx) is a symmetric measure on R (v(dx) is the weak limit of 
E [/iSjv], which is symmetric since the spectrum of is symmetric with respect to 
almost surely), we see that: 

lim -j-E [ti(z 2 I N - AVAOv)" 1 ] =-( ^ (4.103) 

AT^-oo N 1 J Z J R Z - X 

(4.104) 

z 1 — x 
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This implies that, for each zeC\l, 



lim 1e [tr(zI N - X^)" 1 ] = / V ° {X) ~ {dx) . (4.105) 

AT-kx>jV J R Z-X 

and thus, the probability measure E[//^ w ] converges weakly to the measure v o 
i. /•-') \dx). 

Proof of theorem 3.2 ii): using relation (2.12) and lemma 4.9, it is plain to check 
that L(z — x)~ l jiB N {dx) converges in probability to the Stieltjes transform of the 
probability measure 2qj (1 + q)v(dx) + (1 — g)/(l + q)S (dx). This convergence holds 
for finite dimensional vectors (L{zi — x)~ 1 fj,B N {dx)), i = 1, ■ ■ ■ ,d) as well. Using the 
fact that the set of functions {(z — x)~ x ,z 6 C \ M} is dense in the set C*o(M) of 
continuous functions on K. going to at infinity, we can show, for each / 6 C (IR), 
that J f(x)fiB N (dx) converges in probability to J f(x)(2q/(l + q)v(dx) + (l — q)/(l + 
q)6 {dx)). But, since h Bn ( r ) = 2 <?/(l + Q) V ( R ) + i 1 ~ <?)/(! + <?)<$o(K) = 1, this 
vague convergence can be strengthened in a weak convergence. With the relations 
H B 2 n = 2N/(N + T)h Rn + (T - N)/{T + N)5 and the fact that / f{x)n B 2 N {dx) = 
J f(x 2 )fiB N (dx), it is plain to conclude that (jlr n converges weakly in probability to 
v o (x 2 )~ 1 (dx). 

Proof of theorem 3.2 Hi): again using relation (2.12) and lemma 4.9 together with 
Borel-Cantelli's lemma, one can show that the two spectral measures ^b n converges 
weakly almost surely to 2qj (1 + q)v(dx) + (1 — g)/(l + q)5 (dx). It is then easy to 
deduce as before that [1r n converges weakly almost surely to v o (x 2 )~ l (dx). 

A. Auxiliary lemmas 

Lemma A.l. Let A be a n x n complex matrix such that the Hermitian matrix 
M = AA T has spectral radius X ma x- Then, for all i, we have: 

n 
3=1 

Proof. It is straightforward to see that all the entries of M are, in modulus, smaller 
than \ max . On the other hand, we have: 

n 

M ll = J2\ A i3 I 2 • 

3=1 

and, thus: 

n 

^ I A-ij | 2 si A max . (A. 2) 

3=1 

□ 
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Lemma A. 2. There exists C > such that for each N e N and k e {1, . . . , iV}: 



E 



iV 



^^( S )r,(t)Gr)(z) si 



C 



Similarly, for each iVeN and fee {1, . . . , iV} ; i e {1, . . . , N}, we have the following 
inequality concerning the conditional expectation with respect to M l : 



E 



N 



J2 n(sMt)G^\z) st 



\M l 



C 



Proof. We first expand the square and use the independence of (r k (s)) s from G$ +k \z): 



E 



N 



8±t 



TV 



Now we compute 

E[r k (s) 2 r k (t) 2 ] =E 



2j2^[n(s) 2 r k (t) 2 ]E 



M fc (^^, —)M k < 



G 



(JV+fc) 
TV 1 



/V 'TV 



max ( 1 



s -1 / i-1 

N N 

J_ J. 

JV /" JV 



T 



V>(2) 



max ( 1, 



\r — u\ 

T 

\r — u\ 



drdu 



</>(2) 



drdu 



We consider N large enough so as to make 2/N ^ r. The above integral is then 
plain to compute and we get 



y(2)(o2-V>(2) _ 2 n 
E [r fc (s)W] ( 



^(2))(2-^(2))A^ 2 -^ 2 ) 
Thus we have for some positive constant C 



(A.3) 



E 



A? 



^r fe ( S )r fc (t)G^ +fc) (^) st 



C 



N 



AT2-V(2) 



C 1 



U(iV+fc)/ N 



N i-m |9f(^)|2' 

where we have used the fact that almost surely: 

27V 

1 

2iV~ 



1 27V 



s,t^N+k 



2 1 
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It just remains to see that '0(2) = 7 2 . To prove the second relation, we follow the 

r N 



same argument by noticing that (ri(t))t and G^ ,N+L \z) are independent condition- 



ally to M l 



□ 



Lemma A. 3. There exists some constant c > such that for each N G N and 
/■ ■: (1 V}.- 

N 

^ n(A;)rj (A;)G^ } (2) n+ln+j 

i+3 



E 



TV" 



Proof. Again we expand the square and we use the fact that, conditionally to the 
(M ! )j, the quantities ri(k),rj(k), G^\z) N+i ^ N+ j are independent and ri(k),rj(k) are 
centered. Indeed, conditionally to the (M ! )„ the variables ri(k),rj(k), G^\z)N+i,N+j 
involve different increments of the Brownian motion. Thus we have 



E 



N 



( z )N+i,N+j 



N 

= J2^Hk) 2 r J (k) 2 ]E 

N 

< j2nn(k) 2 ]E[ rj (k) 2 ]E 

N 

c 



z ) N+i,N+j 



( z )N+i,N+j 



G n\ z ) N+i,N+j 



where we have used the fact that almost surely: 

27V 

1 

2N 



_ 27V 

— V \g 

T- 1 ^ I 



{k) (z)- ■ 

TV \ z )l,3 



2 1 



|9(*)| S 



Proof of Lemma J^.13. We define the function f^ € on the interval [0, 1] by 

f k /{x) = NM k >\I l N ) ifxeik. 

Notice the relation: 

TV 1 

J2M^I t N )E[G N (z) tt ]= / f k /(r)dE[L]f](dr). 
t=i J o 



□ 
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Then, by stationarity, we have: 

■l 

E ' 1 



f% e (r) dE[L]f](dr) - / e**M dE[L^]{dr) 
o Jo 

N 

USf(r) -(?*<■'>) dE[L)f](dr 



iV 



sup E 



AT 



|9(z)| r 



sup / E 



p/ 1 J/ 1 



1/2 



du 



N 



sup 



(2e^(°)-2e^( r - M )) 1/2 du. 



Because of the continuity of the function p e over [0, 1], we have 



E 



1 f k /(r) dE[L]f](dr) - f e**M dE[L^](dr) 
o ./o 



In a quite similar way, we can prove that 
E 



and 



E 



e^*<f) p {r) dE[L]f](dr)- / e^ (r) dE[L^]{dr) 
Jo 



e"**<f) p (r)K z (r)dr- I e^ {r) K z {r)dr 



as iV -> oo. (A.4) 



—7-0 as p — y oo uniformly w.r.t. N 

(A.5) 



—7-0 as p — y oo uniformly w.r.t. N 

(A.6) 

where (</> p ) pe N is a regularizing sequence and * stands for the convolution. Fur- 
thermore, for each fixed p and because of the weak convergence of E[L^ 2 ] towards 
K z (x)dx, we have almost surely 

/ * <f) p {r) dE[L)f](dr) / * <f) p {r)K z {r) dr as N ^ oo. (A.7) 
Jo Jo 

We prove the result by gathering (A.4) (A.5) (A.6) and (A.7). □ 



B. Sup of MRW 

Here we prove 

Proposition B.l. We have for all k — 1, . . . , N + 1 

(In AO 2 



E 



sup r k (ty 

t=l,...,N 



< c 



C(2a)-1 



for some positive constant C. 
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Proof. To prove the result, we first prove 

Lemma B.2. There exists a constant C such that, if (Xi)i ^ , <g n are iid centered 
Gaussian random variables then: 



E 



max IXA 4 



€ C 



max E[X 2 l 2 (lniV) 2 . 



Proof. By homogeneity, it suffices to assume that E[X 2 ] = 1. Then we have for all 
5>0 



E 



POO 

max \XA 4 + N FilXA 4 > t)dt 

POO 

^5 + 2N J F(Xt > t 1/4 )dt 



€ 5 



€5 + 



2N 



'2ix Js 

an r 



'2ix Jy/S 

AN 



-^dt 



e'Hdt 



and this last expression can be made smaller than C(lniV) 2 by choosing 5 = (IniV) 2 . 

□ 

We want apply the above lemma after conditioning with respect to the law of 
the MRM M k : 



E 



sup r k (ty 

t=l,...,N 



E 



E 



sup r k (t) 4 \M k 

t=l,...,N 



Notice then that, conditionally to M k (0, —) = x%, . . . , M k (^^-, 1) = x^, the vector 
(rfe(l), . . . , r k (N) has the same law as the increments of B: (B X1 — B Q , . . . , B XN — 
B XN _ 1 ). By applying Lemma B.2, we deduce that 



E 



sup r k (t) 4 \M k 

t=l,...,N 



Thus we deduce 
E 



sup r k (ty 

t=l,...,N 



^C(lniV) 2 max M k (- -, — ) 2 . 

V ; t=i,...,N ^ N N' 



max M k ( t — -, 4) 
t=i n v N N' 



<: C(lniV) 2 E 

Finally we have for all 5 > and for a > 1 such that ({2a) > 1: 
E 



(B.l) 



max M fe (— ,-) 
=i,..,n v N ' N J 



pOO r 1 ± 

^5 + N F(M k (— ) 2 >x )dx 



^ 5+N ri E [M k (y4) 2 i dx 

L x a 1 V N ' N> J 



N ' N' 
- 1 

~N~' N- 



^5 + C5 l - a N 1 ~ a2a) 
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for some constant C only depending on a, r and r ) 2 . Choose now 5 = N ^ ' so as 
to get 



E 



sup r k (t) 4 



< (1 | C ) (lniV)2 

iV a 



(B.2) 
□ 



C. Girsanov transform 

Lemma C.l. Let /x fre an independently scattered infinitely divisible random measure 
associated to (if),0), where 

Wq e R, ip(q) = mq + \o 2 q 2 + / (e qz - l)v(dz), 

< +00 and ^(l) = 0. Let B be a bounded Borelian set. We define a new 
probability measure ¥ B (with expectation K B ) by: 

VAmeasurable set, F B (A) = E[l A e^ B) }. 

Then, under F B , ji has the same law as fi+fi B where ji B is an independently scattered 
infinitely divisible random measures independent of /1 and is associated to (ip B ,9 B ) 
given by 

fo(q) = qa 2 + I (e qx - l)(e x - l)v{dx) 
Jr 

9 B (-) = 9(-nB). 

Proof. It suffices to compute the joint distribution of p disjoint sets A\, . . . , A p . We 
have for any Ai, . . . , A p e K: 



E, 



3 Ai/i(Ai)+-+Ap/i(A p ) 



= E 
= E 
= E 
= E 



Xifi(A 1 )+-+X pl j.(Ap)+fi(B) 



^Xi l x(A 1 \B)+-+X P ti(A p \B)+X 1 [i(A 1 nB)+-+Xp tJ ,(A p nB)+n(B) 



X 1 fi(A 1 \B)+---+X p ^A p \B)+(X 1 + l)fi(A 1 nB)+---+(X p +l) f i(A p nB)+fi(B\\J , l =1 A i ) 



: ,X lt i(A 1 \B)+-+X pl i(A p \B) 



E 



3 (Ai+l)/i(AinB)+-+(A p +l)/i(A 



P nB)j 



e i>(X 1 )8(A 1 \B)+---+i;(X p )e(A p \B) e i;(X 1 +l)e(A 1 nB)+---+^(X p +l)e(A p nB) 
e ^(Ai)e(Ai)+-+t/>(Ap)e(A p ) e (V(Ai+l)-^(Ai))e(AinB)+-+(V(A p +l)-V(A p ))e(A p nB) 



Then it suffices to notice that: 



i[>(q + I ) - i[>(q) = m + a 2 q + \a 2 + / I 

2 Jr 



e 9 *)z/(dz) 



and V(l) = 0. 



□ 
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Lemma C.2. If the process u € is defined as u e (x) = fj,(A e (x)) where \l is an indepen- 
dently scattered random measure associated to (ip,6) with tp(q) = —iq r y 2 /2 — q 2 ^ 2 /2 
and 9 given by 2.3, then: 



limE 



u) e (x) 



r)e UJ ^ r ^dr 



E 



r — x\ 



-i 



K z (r)M{dr) 



where M is the lognormal MRM. 



Proof. One can check that (u € (x)) xe [o ; i] is a stationary gaussian process with covari- 
ance given by 7 2 p e (x — y). So, using Girsanov transform, we can write: 



E 



z - C K z (r)e u *( r )dr 



E 



K z (r)e l2Mr - x) e^ {r) dr 



We are interested in the limit when e goes to of this latter term, we thus 
approximate it with a simpler term: 



E 



z- # z (r)e 72 * (r - x W r > dr 



E 



o 



: - / K z [r) ( -^—Y e^ {r) dr 
to \\r - x\ 



-i 



E 



\K z (r)\e^ {r) 



3 7 2 p £ (r-s) _ 



3 7 2 p E (r-x) _ 



r — x\ 



dr 



\r — x\ 



dr 



(C.l) 



where we have used Lemmas 4.4 and 4.8 and the normalization = 0. 

Because r y 2 < 1, the dominated convergence theorem implies that C.l converges 
to when e goes to 0. 

We thus look at the limit when e goes to of the term: 



E 



to \\r - x\ 



3 w ' (r) dr 



The random variable 



K z (r) 



r — x\ 



M(dr) 
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is well defined and is finite almost surely since: 



E 



KJr) 



o 



T 



7 



\r — x\ 

And thus, we can compute: 
• i 



M{dr) 



\K z {r) 



T 



r — x\ 



dr < +oo. 



E 



z — 



KJr) 



o 



r — x\ 



E 



z- \ KJr) 



\r — x\ 



M{dr) 



E 



K z (r) 



\r — x\ 



>- {r) dr - M{dr)) 



and, for all n6N, this latter term is smaller than 



E 



+ E 



+ E 



KJr] 



T 



\r — x\ 



mm 



r — x\ 



n 



' {r) d7 



K 7 (r) min 



KJr) 



T 







\r — x\ 

,2 



n) (e We(p) dr- M(dr)) 



\r — x\ 



7 



mm 



r — x\ 



. n 



M{dr) 



The two quantities C.2 and C.4 are smaller than 

. 2 



\Kjr) 



o 



r 



r — x\ 



mm 



T 



\r — x\ 



dr 



(C.2) 
(C.3) 
(C.4) 



(C.5) 



and thus converge to 0, uniformly in e as n goes to infinity. 

2 

For a fixed n, the function min((r/|r — x|)+ , n) is measurable and bounded and 
thus it is plain to see that, for a fixed n, the term C.3 goes to when e goes to 0. 
The lemma follows gathering the above estimates. □ 

Lemma C.3. If the process u 6 is defined as u t (x) = fi(A 6 (x)) where \i is an indepen- 
dently scattered random measure associated to (<pJJ) where <p is given by (2.1), i.e. 



ip(q) = imq 
and where 9 given by (2.3), then: 



limE 

e->0 



3 OJ e (x) 



z - Jg 1 K z {r)e^ r )dr 



Y q + 



E 



iqx 



l)v{dx) 



lr — x\ 



-i' 



K z (r)Q(dr) 
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with k = 7 2 + j R {e x — l) 2 v{dx) and where the random Radon measure Q is defined, 
conditionally on a MRM denoted by M whose structure exponent is ((q) := q — 
(p(—iq), as the almost sure weak limit as e goes to of the family of random measures 
Q e (dt) := e Wt< $ 'M(dt) where, for each e > 0, the random process u e is independent 
of M and defined as u t (t) = ~fi(A t (t)) where ~p is the independently scattered log 
infinitely divisible random measure associated to ((p,9(- n A (x))) where 

<p(p) = iptf -k)+ f (e ipx - l)(e x - l)v(dx). (C.6) 
Jr 

Proof. We want to apply Lemma C.I to the process u e . If we set B = A e (x), Lemma 
C.l tells us that, under P_g, the process u e possesses the same law as the process 

w W( r ) + w ( 2 )( r ) W ith w e (1) (r) = fi {1) (A e (r)) and uf\r) = /i (2) (A e (r)), 

where independent independently scattered log infinitely divisible ran- 

dom measures respectively associated to (</?, 9) and ((p^ 2 ',9^ 2 ') with: 

^){ q ) = %1 2 q + ! ( e iqx - l)(e x - l)v{dx) and 9 {2 \-) = 9{- fl A e {x)). (C.7) 
Jr 

Define: 



k = 7 2 + / (e x - l)' 2 v(dx),p(q) = ¥ {2 \q) -iqK,ip(q) = <p(-iq). 



(C.l 



Notice that ip is then normalized so as to make ip{l) = V'(O) = 0. Let us define the 
process uJ e by: 

uj e ( r ) = ujf\r) - K9(A e (r) n A e (x)) = J 2 \r) - np e (r - x), (C.9) 

and notice that E[e iqJJ ^] = e ^ q ^ r - x \ 
We can now apply Lemma C.l: 



E 









= E 


(- 


/ K z (r)e^ {r)+Kp ^ r 






Jo 



We are interested in the limit when e goes to of this latter term, we thus 
approximate it with a simpler term: 

x -r 

e U E (r)+oj E (r)+Kp E (r-ic)^ ^ r 







E 


N 







E 



\%(z)\ 2 



E 



uj e (r)+u> e (r) 



\r — x\ 



K z (r) dr 



l 



3 Kp E (r-x) 



r — x\ 



\K,(r) \ dr 



\r — x\ 



dr 



(CIO) 
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where we have used Lemmas 4.4 and 4.8, the normalizations "0(1) = 0, ^(l) = and 
the independence between UJ e and u t . 

Let us show that k < 1. Indeed, we have: 



7 2 + / (e x - l) 2 u{dx) 
Jr 

7 2 + / (e 2x - l)v{dx) - 2 / (e x - 



V 2 + / (e 2x -l)u(dx) + 2(m + ]p 2 ) 



= 2m + 2 7 2 + / (e 2 * - l)u(dx) 
Jr. 

= m 

where, in the third line, we used the fact that ^(l) =0 (which implies the relation 
j R {e x — l)u(dx) = — (m + 7 2 /2)). We will now show that '0(2) is strictly less than 
1. It suffices to show that £(2) > 1. Using the concavity of the function (, we have 
the inequality: 

c(2 + t )-cd) <c(2) _ c(1) (au) 

1 + e 

and with assumption 3.7, we see that ((2) — ((1) = ((2) — 1 > 0. We can thus 
conclude that k < 1. 

Because k < 1, the dominated convergence theorem implies that CIO converges 
to when e goes to 0. 

For each Borelian set A of [0; 1], the family M e (A) := f A e UE ^dr,e > is a 
positive martingale with respect to e and that it converges almost surely to M(A). 
With the assumption 3.7 and in particular the condition ((2 + e) > 1, we can show 
(see [1] for a proof) that the family (M t (A)) t>0 is in fact uniformly integrable. In 
particular, if we let T t be the sigma field generated by the family of random variables 
(a; r? (r)) r)>eire iR, we have the following almost sure equality: 

E[M(A)|J" e ] = M e (A). (C.12) 

Conditionally to the random measure M, the family P £ (A) := j A e u ^M(dr), e > 
is also a positive martingale with respect to e. Thus, P e {A) converges almost surely 
to a random variable that we will denote by P(A). We know that this defines a 
random Radon measure P on [0; 1] and that the family of random Radon measures 
P e converges, when e goes to 0, weakly almost surely to P in the space of Radon 
measures. Denote, conditionally to the random measure M, by Pm the law P[-|M] 
and let us show that the family (P e ([0; l])) e >o is P^-uniformly integrable. Let 6 be 
such that -0(1 + 5) < +oo (we can show, using the condition ip(2 + 5) < +oo, that 
that there exists such 5 ). We will show that the family (P e ([0; l])) e >o is uniformly 
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bounded in L 1+s (Fm)- Indeed, conditionally to the random measure M: 



E 



] M(dr 









> 




11 









;(i+*)«.«JVf(d r ) M[0; if 
:^(i+3)p e ('-*)M (dr)M[0; if 



^ M[0; lfe 



SMl+S) 



r — x 



M(dr) < +00. 



The family (-P e ([0; l])) e >o is therefore P A /-uniformly integrable, in particular, P e ([0; 1]) 
converges to P([0; 1]) also in L 1 , which implies that P is a non degenerated random 
measure. Moreover, denoting by J- \ the sigma field generated by the family of ran- 
dom variables (cJ r? (r)) r)>ergK , we have, almost surely, conditionally to M, for all 
Borelian set A of [0; 1]: 

E M [P(A)\F e ] = P e (A). 

Now, as before, it is easy to see that the family Q e (A) := J A e^OH^O)^ e > 
is also a positive martingale with respect to e. Therefore, Q e {A) converges almost 
surely to a random variable that we will denote by Q{A). This defines a random 
Radon measure Q and the family of random Radon measure Q e converges, as e — > 0, 
weakly almost surely to Q in the space of Radon measure. We want to show that 
the two random measures P and Q have the same law. 

Gathering the above arguments, we can write, almost surely: 

E [P(A)\a(F e ,? e j\ = E [E[P(A)\T e ]} 

E I e^ (r) M(dr)\F t 



and the latter quantity has the same law as Q e (A). Since the martingale (E[P(A)\a(J : ' e , ^ r e )]) e>0 
is uniformly integrable, we deduce that the family {Q t {A)) e>0 is also uniformly in- 
tegrable. Hence, both random variables P(A) and Q(A) have the same law. We can 
show easily that in fact the two random measures P and Q have the same law. In 
particular, Q is non degenerated. 

It is now easy to see that, for all bounded and continuous function /, the two 
random variables J R f(r)P(dr) and L f(r)Q(dr) have the same law. By regularizing 

the function ( -r 2 — - ) and with the dominated convergence theorem, we conclude as 



|i — x\ 

in the proof of lemma C.2 using the fact that k < 1 that: 



Kjr) 



Q(dr] 



(law) 



K z {r) 



x\ 



P{dr). 



(C.13) 



Gathering the above argument and letting e go to concludes the proof. 



□ 
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